LU TP 98-14 
jhep-ph/9807541 
July 1998 



Is there screwiness at the end 
of the QCD cascades? 



Bo Andersson, Gosta Gustafson, Jari Hakkinen, 
Markus Ringner and Peter SuttonF] 

Department of Theoretical Physics, Lund University, 
Solvegatan 14A, S-223 62 Lund, Sweden 



We discuss what happens at the end of the QCD cascades. We show that, with just a few 
reasonable assumptions, the emission of soft gluons is constrained to produce an ordered 
field in the form of a helix. We describe how to modify the Lund fragmentation scheme 
in order to fragment such a field. Our modified fragmentation scheme yields results which 
are consistent with current experimental measurements, but predicts at least one signature 
which should be observable. 



1 bo@thep.lu.se, gosta@thep.lu.se, jari@thep.lu.se, markus@thep.lu.se, peter@thep.lu.se 



1 Introduction 

In QCD the production of two colour charges which subsequently move apart will lead to the 
production of further colour radiation. This can be described in terms of the fundamental 
field quanta, the gluons, but it is also possible to describe the ensuing radiation in terms of 
dipoles. This property arises because in non-abelian theories the emission of an extra gluon 
from a gluon-gluon dipole can (to a very good approximation) be modelled as the destruction 
of the original dipole and the creation of two new dipoles. In this way the change in the 
colour field can be described as an increasing cascade of dipoles. The end of this cascade 
occurs when the dipole masses are so small that helicity conservation prevents further real 
gluon emission. In this paper we examine what happens at the end of this cascade. We find 
that the conditions are favourable for the field to utilize the azimuthal degree of freedom 
and wind itself into the form of a helix. This corresponds to a close-packed configuration of 
gluons in rapidity-azimuthal-angle space. 

We begin by describing a toy model which contains the relevant features, namely a ten- 
dency to emit as many gluons as possible and the constraint that gluons are not too "close" 
to each other (which arises from helicity conservation). In this simple model it is clear that 
at the end of the cascade an ordered field emerges with the characteristics of a helix. To 
progress beyond this model we use the Lund model of QCD. In the Lund picture hard glu- 
ons are represented as excitations of a relativistic string which connects a quark, anti-quark 
pair. However, the gluons from which the helix is built up are too soft to be modelled in this 
way. Instead we introduce a helical semi-classical field and thus develop a modifed version 
of the Lund fragmentation scheme. Our modified fragmentation scheme enables us to study 
whether the consequences of a screwy field can be detected in the final state particles. We 
find that if events with hard gluons are excluded then the screwiness of the field may be 
observed. 

2 The dipole cascades; increase and decrease of phase 
space 

In order to describe what can happen at the end of the QCD cascades we will provide a 
brief description of the cascades. We will in particular discuss the consequences of helicity 
conservation in the emission of partons. 

The well-known formula for dipole emission of bremsstrahlung is 



where a is the effective coupling, k±, y, and are the transverse momentum, rapidity and 
azimuthal angle respectively, although the azimuthal angle dependence is usually neglected. 
The final factor, corresponds to the spin couplings. We will briefly consider the precise 
definitions before we consider the implications. The effective coupling for QCD in the case 
of a gluonic dipole is given by 
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The occurrence of the number of colours, N c , and the factor 1/2 in the QCD coupling is due 
to early conventions, whereas the result that the running is governed byl/c = 6/llisa basic 
gauge group independent result. It only depends upon the fact that in non-abelian gauge 
theories there is a three-particle coupling between vector particles, e.g. the colour-8 gluons 
in QCD. (The four-gluon coupling also occurs to preserve the symmetry, but it does not play 
a role in this connection). We neglect the flavour term — 2n//3 which should accompany 11 
in the denominator because it is a small effect related to the possibility of gluon splitting; 
g -> qq. 

The transverse momentum and the (dipole cms) rapidity are defined in a Lorentz invariant 
way in terms of the squared masses of the final state partons (the emitters are conventionally 
indexed 1 and 3 and the emitted field quantum 2): 

Sij = (h + kj) 2 = 2kikj = 2k_ak±j [cosh (Ay) j.,- — cos(A0)jj] 

S = S 12 + S 2 3 + «31 

S12 = s(l ~X 3 ), S 2 3 = s(l -Xi) 

,2 S12S23 

y = hogM (3) 
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Here X\ and £3 are the final state cms energy fractions of the emitters. Requiring energy 
momentum conservation limits the allowed emittance region to 

fc ± cosh(y)<^. (4) 

This region can conveniently be approximated as \y\ < (L — k)/2 with the variables L = 
log(s/A 2 ) and k = log(Ar[/A 2 ). This means that the (approximate) phase space available 
for dipole emission is the interior of a triangular region in the (y, /«)-plane with the height 
and the baselength both equal to L. The inclusive density inside the triangle is, in this 
Leading-Log Approximation (LLA), given by the effective coupling a according to Eq. ([]]). 
The rapidity range, L — n, is of course the length of a hyperbola spanned between the emitters 
in space-time (or energy-momentum at the scale kj_). 

If we consider an initial qq dipole emitting a gluon then the probability for the produced 
qgq system to emit a second gluon is a complicated expression QXJ] . In case the transverse 
momenta of the first and second gluon are strongly ordered, k±i 3> k± 2 , it is a very good 
approximation to treat the second emission as independent emission from two dipoles 0. 
For an exclusive statement, for example the probability dP to emit the first gluon with a 
certain it is necessary to multiply the inclusive formula in Eq.(|l]) with a Sudakov 

form factor A s containing the probability not to emit above Ki, 

A a (L,/ci) = exp(- / dn), 
c^P(q,gi,q) = dn(Ki,yi)A 5 (L,«i). (5) 

The probability to emit two gluons is then, in the approximation that the second gluon is 
emitted by two independent dipoles, given by 

<£P(q,gig 2 ,q) = rfP(q,g 1; q) [dP(q,g 2 ,gi) + rfP( gl , g 2 , q)] (6) 
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in easily understood notations. The approximation in Eq. results at most in a percentage 
error over all phase space ||. Thus, contrary to QED where the chargeless photons still 
leave the e + e~-current as the single emitter, the 8-charge gluon (gj in QCD changes the 
original qq dipole into two dipole emitters, one between q and g x and one between g : and q, 
and each can independently emit the second gluon (g 2 ). The requirement for the validity of 
the approximation in Eq.(^) is that k±i > k± 2 or else the indices are exchanged. 

The two independent dipoles are moving apart (with g l as the common parton). This 
means that they have together a larger effective rapidity range for the emission of g 2 , i.e. 
the original hyperbola length L = log(s) is exchanged for two hyperbolas with the combined 
length log(s qgl ) + log(s gi q) = L + \og(k'^ 1 ). From any one of the two new dipoles we may 
then emit the second gluon, thereby producing three independent dipole emitters and the 
process can be continued towards more dipoles; ordering the process in k± downwards. The 
available phase space for further emission is increased after each emission, as can be seen 
from the increased total length, L, after the first emission. This description of the QCD 
cascades is called the Dipole Cascade Model (DCM) 

We will now consider the polarisation sum contribution in Eq.(|l|). Its precise properties 
depend upon whether we are dealing with a qq, qg or a gg dipole, but it stems from the 
spin couplings between the emitter(s) (it is essentially sensitive only to the closest emitter) 
and the new field quantum. These couplings contain the property that helicity is conserved, 
which is true for all gauge theories. This means that if a spin-1/2 parton emits a spin-1 
parton, the spin-1 parton must go apart from the emitting particle in order to conserve 
helicity and angular momentum. They have to go even further apart in the case of a spin-1 
parton emitting a spin-1 parton. To estimate the separation we consider (for fixed k± (or 
k)) the available rapidity range: 

/ ^ dy = L- k- c + 0(k 2 Js) (7) 

2/min 

where c = (11/12 + 11/12), (3/4 + 11/12) or (3/4 + 3/4) depending on whether the emitters 
are gg, qg (gq) or a qq dipole ||. The quantities c are written as sums to show that a spin-1 
(g) emitter and a spin-1/2 (q or q) emitter has an empty region surrounding it in rapidity 
of size 11/6 and 3/2, respectively. In order to obtain this result we note that in terms of the 
x-variables introduced in Eq.(|^) the factor \1> is (x" 1 +x% 3 )/2 with a+3 = 1 — kj_ exp(±y) / y/s 
and 711,3 equal to 2 or 3 for q(q) and g, respectively. ?/{max,min} are determined from the 
energy momentum requirement in Eq.(|4]). 

A note of caution should be issued at this point. For given s and k± there are two definite 
limits in rapidity y min < y < y max , and there is then a depletion of emissions due to helicity 
conservation, in regions close to y min and 2/ max - It is in general a poor approximation to put 
the factor \1/ to a unit stepfunction for y min + c/2 < y < y max — c/2 although it works when 
the rapidities and azimuth are integrated out. A closer examination provides a y-distribution 
with similarities to a finite temperature Fermi distribution. We will nevertheless refer to this 
feature as "the excluded region" around each gluon. 

We note that in the process g — > qq, where the spin-1 parton emits two spin-1/2 partons, 
that the fermion pair "prefer" to be parallel, since there are no poles in this decay distribu- 
tion. However the process g — > qq is suppressed compared to the process g — > gg and is in 
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general neglected. The DCM will in this way produce a fan-like set of dipoles, which in the 
LLA increases the phase space (the total available effective rapidity range) for further emis- 
sions. However, including the influence from the polarisation sum (which is essentially the 
approximation scheme called Modified LLA) there is in each emission also a depleted region 
around an emitted parton, in practice c = 11/6, because the gluons completely dominate the 
process. At large energies, but not too large transverse momenta, one may in general neglect 
the restrictions but they will be very noticable at the end of the cascades. For example, with 
a dipole mass of 3 GeV the typical rapidity range available for gluon emission is about 4 
units, and it is then very noticable to exclude 11/6 units. 

It is interesting that the average region excluded due to helicity conservation also occurs 
in connection with the properties of the running coupling. To be more precise, we consider 
a change of scale in the definition of a field quantum and its interaction. A change of scale 
means that the field operator, which has been normalised to a single quantum at one scale, 
and the coupling constant, which likewise has been normalised at the original scale, will 
both change. These changes can be read out from the Callan-Symanzik equations and the 
/5-function contribution, stemming from the change in the coupling constant, can be written 
as 



where a change in a quantity A4, when the observation scale is decreased from the level 
k = \og(k\) to k — d,K, is considered. The decrease accounts for the minus sign on the left 
hand side. According to the DCM there is then at this new scale not only the possibility to 
emit new gluons but also, at the next order in the coupling a s , the possibility to reabsorb 
already emitted gluons. 

The operator a s d/da s works like a number operator, i.e. for any function M. = J2 ot™m n it 
provides the number n of possible insertions. The quantities N c a s /2ir and nfa s /4n are the 
couplings for gg — > g and qq — > g (and the inverse processes) while 11/6 and 2/3 corresponds 
to the effective (generalised) rapidity ranges available in these reabsorption processes for a 
given k. It should be noted, however, that this interpretation is gauge-dependent; in almost 
all gauges there are contributions to the /3-function from the vertex corrections. However, 
for a particular gauge choice with the propagator given by —{g^ v — <ik^k u /k 2 )/k 2 , the vertex 
contributions vanish. A closer analysis reveals that the major effect stems from the so-called 
Coulomb gluons, i.e. a charged particle like a field quantum in a non-abelian theory is always 
accompanied and interacts with its own Coulomb field. The 11/6 can therefore be considered 
as the region around the gluon containing its accompanying field. This has been utilized 
for an approximation of the QCD cascades where the available phase space for emission is 
discretised |J. 

3 A toy model for the end of the cascades 

After several gluon emissions there are a set of dipoles with small masses, and there are in 
general very many Feynman graphs which may contribute. The largest diagramatic contri- 
bution is chosen according to coherence conditions in the cascade; in the Dipole Model Q by 
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an ordering of the gluon emissions in transverse momentum, and in the Webber-Marchesini 
model 0, and the model implemented in JETSET ||, according to a choice of kinematical 
variables that fascilitates a strong angular ordering of the emitted gluons. Results from the 
cascade models are essentially equivalent, at least as long as sufficiently hard gluon emission 
is considered. 

The ordering of emissions in the models will lead to dipoles with small masses emitting 
softer gluons. These soft gluons have a transvere momentum, k±, of the same order as their 
emitter and recoils play an important role. At present there exists only a minor knowledge 
of how the recoils should be distributed among the emitters. A sufficiently large recoil on 
one of the emitting (soft) gluons will in general imply that the chosen order is no longer 
in accordance with the coherence conditions. Emitting soft gluons will evidently lead to a 
situation where several, or even very many, paths to the final state are important, and many 
different Feynman graphs may contribute and interfere. 

To investigate the emission of soft gluons we propose a toy model with the following two 
properties: 

I We assume that the effective coupling a is large enough so that there is a tendency to 
emit as many gluons as possible, essentially with the same k±. 

II We assume that the emissions fulfil the requirement of helicity conservation; this implies 
that two colour-connected gluons cannot be closer than a "distance" d = c. 

We will use the following combination as the probability for a given colour-connected 
multi-gluon state 

P=U— ( 9 ) 

1 s jj+i 

where Sjj+i is the dipole mass between the colour-connected gluons j and j + The factor a 
corresponds to the product of the coupling and the relevant phase space region, and (3 to the 
restrictions from helicity conservation, i.e. the requirement of a suitable distance between 
the emitted gluons. Neglecting recoils, we obtain for any order of the emissions in the DCM, 
that the product of factors l/sjj + i can be written in terms of the invariant dipole transverse 
momenta as 

(10) 

where k±j denotes the invariant k± of the dipole from which gluon j is emitted. Eq.@ is 
therefore a simple generalisation of Eq. ([!]). 

The dipole mass can be written as 

Sjj+i = £r[2[cosh(A?/) — cos(A0)] 

~ fci(Ay 2 + A0 2 )[l + (Ay 2 -A0 2 )/12]. (11) 

For simplicity we have set the transverse momenta of the two gluons to be identical. Ay and 
A(f> are the differences between the colour- connected gluons in rapidity and azimuthal angle, 
respectively. We are now in a position to define precisely what we mean by "distance" . We 
therefore introduce a distance measure, d, which is related to the dipole mass by 

d jj+i = \J s 3,j+i/ k l ■ ( 12 ) 
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When the dipole mass and the rapidity region are large the azimuthal dependence can be 
neglected and d ~ Ay. 

The emission of soft gluons has thus been reduced to the following problem; given a certain 
rapidity range and the full accompanying azimuthal range < <fi < 2n how are the colour- 
connected gluons distributed in phase space in order to obtain a maximum of P in Eq.(^), 
keeping in mind that the gluons cannot be too close? 

From Eq.(^) we see that the magnitude of a controls the relative probability between 
different gluon number states. If a is sufficiently large the number of emitted gluons will fill 
the available phase space, and P becomes maximal when the gluons align along a straight 
line in phase space. This helix-like structure is the optimal configuration irrespective of the 
size of a, or of the number of emitted gluons. For a given multi-gluon state there are many 
possible ways to colour-connect the state, where the helix is only one of the possibilities. It 
is of course possible that the sub-optimal configurations are the important ones and swamp 
the helix-like contribution, but there are also many contributions close to a perfect helix. 

We have carried out a numerical study to test whether the contributions from helix-like 
structures survive the phase space effects. Our program calculates all possible configurations 
on a discretized (y, 0) phase space taking into account that gluons must not be closer than 
c to each other. The number of possible configurations grows factorially with the number 
of gluons, but the number of gluons is restricted by the available phase space. We have 
studied a reasonable phase space size of three units of rapidity using a closest gluon-to-gluon 
distance c = 11/6 in all the calculations. Since the fluctuations in dipole k± are limited 
within a narrow range at the end of the cascades and the dependence on dipole k± in Eq. (TJ) 
is rather weak, we set the transverse momenta of the gluons to be constant. 

In Fig.|l| we show the most probable five and six gluon configurations. The points cor- 



responding to a given mass correspond to ellipse-like shapes (Jcosh(5y) — cos(50) ) and in 
order to minimize the distance between adjacent gluon emissions these ellipses must be dis- 
placed so that they correspond to a helix-shaped configuration. The case shown corresponds 
to the optimal situation where it is favourable to "close pack" the gluons irrespective of the 
size of a. 

Taking into account all possible configurations we obtain a distribution in D 2 = djj+i 
which is very broad, cf. Fig.^], but weighting each D 2 with the corresponding P from 
Eq.(^|) we obtain a large and narrow peak close to the most probable colour-connected 
configuration indicating that the gluon configurations have short strings close to the optimal 
helix structure. 

Now that we have established that short strings are preferred we investigate in more detail 
if they are helix-like in general. To this aim we will introduce a new possible observable, 
"screwiness" . At this point it is only a theoretical observable, but later on we will show how 
to use it for the final state hadrons. We define screwiness S from the values of {yj,<f)j) for 
the emitted gluons in accordance with the toy model, 



The first sum is over all the configurations e found in the phase space and the second goes 
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(13) 
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Figure 1: The most probable configurations with five and six gluons using c = 11/6. (The 
cylindrical phase space has been mapped onto a plane). The gluon exclusion region for each gluon 
is indicated with the ellipse-like shapes. The line segments show the colour field and should form 
a straight line for a perfect helix. The discrepancy is due to the discrete phase space used in our 
numerical analysis. 

f(D 3 ) F 1 



0.15 - 



0.10 - 



0.05 - 




10 20 30 40 50 Tj 2 

Figure 2: The unweighted (dashed line) and weighted (solid line) squared length distributions, /, 
of configurations with six gluons. 

over the gluons in the configuration. For cj-values close to zero, screwiness must be small if 
the gluons are emitted isotropically in the azimuthal angle. For large values of u the phases 
should be close to chaotic and then screwiness only depends on the mean number of emitted 
gluons. 

In Fig.||| we show the screwiness distribution including contributions from all configurations 
with a specific number of gluons. Two cases are shown, firstly configurations with the 
maximum possible number of gluons (in a three unit rapidity phase space this is six gluons), 
and secondly those corresponding to five gluon states (the contributions corresponding to 
even smaller number of gluons show similar distributions). There are two noticable broad 
peaks with their mean values close to u — ±2tt / c. Since the helix structure has no preferred 
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Figure 3: Screwiness in the toy model for five (dashed line) and six (solid line) gluon states in a 
rapidity region of three units with c=ll/6. 

rotational direction the distributions should be even. The small apparent asymmetry is due 
to numerical effects. We have also analysed the configurations for c = 1.5 and 3 and these 
results are independent of the minimum gluon-to-gluon distance. 

From this toy model we see that if we fill the phase space with soft gluons, which are 
forbidden to be too close to one another, then they tend to line up along a helix structure, 
since the colour-connection between the gluons prefer to be as short as possible. 




4 Modelling the helix as an excited string 

In order to consider the consequences of the helix-like colour field which we obtained in 
Section |3| it is necessary to provide observables in terms of the final state hadrons. A first 
attempt to model such a field is to approximate it by the emission of a set of colour-connected 
gluons with the same transverse momentum k±. We may then consider the properties of 
the final state hadrons, as produced by the Lund string fragmentation model. We very 
quickly find that in the competition between increasing the multiplicity versus increasing 
the transverse momentum of the hadrons the model uses the first possibility only. In this 
section we will be content with giving the basic argument for why the helix cannot be 
described as gluonic excitations on the string field. 

Suppose that a gluon with transverse momentum k± is moving transversely to the constant 
(k) force field, then it is possible for the gluon to drag out the string field the distance 
I = k±/2n (a gluon experiences twice the force acting on a quark). On the other hand, in a 
quantum mechanical setting such a gluon is only isolated from the field if the wave-length 
of the gluon A ~ 2ir/ k±_ is smaller than I and therefore 

k\ > 4nK 

f ■ = 2L (u) 
(this is similar to the Landau- Pomeranchuk formation time arguments). From the first line 
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in Eq.(14) we obtain the requirement that a "real" gluon must have a transverse momentum 
larger than k± = 1.6 GeV. 

We conclude that the helix field cannot be described in terms of a finite number of gluon 
excitations on the Lund string. The many small-fcj_ excitations in the model tend to increase 
the final state particle multiplicity (with small fluctuations) rather than to produce transverse 
momentum for the particles. The interested reader can find a more thorough investigation 
of the problems associated with the fragmentation of soft gluons in appendix |A]. 



5 A semi-classical field at the end of the cascades 

We will now consider the possibility that a (semi-) classical colour field is produced at the 
end of the perturbative QCD cascades that cannot be described solely in terms of gluonic 
excitations on the Lund Model string field. The properties of this field should be in accor- 
dance with the toy model that was described in Section [|. Thus the internal colour quantum 
number should be correlated to the external space-time (energy-momentum space) behaviour 
so that the colour field has a helix structure, i.e. the colour field lines are turning around a 
spacelike direction, from now on called the 1-axis. 

We may describe the expected field in terms of a wave-packet of energy-momentum space 
four- vectors, ke, corresponding to the colour current (the index 8 stands for the parameters 
describing the wave-packet). We will assume that the vectors ke always have a constant 
virtuality kg = —m 2 . We further assume that the helix colour field is itself emitted from the 
current as a continuous stream of gluons dk, colour-connected along each emission vector, 
kg. They should be obtained by differentiating the vector ke (we are generalising the physics 
picture from a ladder- diagram as in Fig.f|, where the "propagator" vectors {k}j are emitting 
the gluons dkj = kj — fc 3 -_i). 

The most general description of such a vector is (we use lightcone coordinates along the 
01-direction and transverse coordinates in the 23-plane and we do not worry about the initial 
values) : 

k 9 = rn[cos(8)(exp(y), — exp(— y), 0, 0) + sin(0)(O, 0, cos(cr0), sin(<70))] 

= m[cos(0)ei(y) + sin(0)e_j_i(cr0)]. (15) 

Here m is a constant parameter, y is the rapidity and <fi the azimuthal angle. We have 
introduced a as a constant describing the relative motion in rapidity and azimuth. We will 
put a = 1/2 later in order to get (ft as the azimuthal angle of dk. Finally 9 is the variable 
describing, on the one hand, the size of the fluctuations in the longitudinal and transverse 
parts, and on the other hand, the properties of the wave-packet. 

Assuming that the emitted field quanta dk are massless, we get, 
dk 2 = 

small 6 1^-1 +6 2 =1^1 . (16) 
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Figure 4: A current with constant virtuality, k 2 = —m 2 , emitting massless field quantas, dk 2 = 0. 



We have used the differential dl = ydy 2 + d(a(j)) 2 . Therefore, the assumptions of constant 
virtuality of k and the masslessness of dk imply that the variable 6 should fulfil the pendulum 
equation according to the first line of Eq. ([l6|) . In the limit of small |#|-values this becomes a 
harmonic oscillator equation, assuming that the quantity dyjdl is a (small) constant along 
each vector kg. For consistency we will then make the change di — *> ad<p. This is the second 
line of Eq.(|T6D and using the notation dy/d(f) = r we obtain as a classical description (again 
neglecting the boundary values): 

= -cos(a<j)) (17) 
a 

If we choose a = 1/2 to make <fi the azimuthal angle of dk, then we find that the field 
emission vectors dk/dcj) and the corresponding current vector ke are (in the approximation 
of small oscillations): 

dk 

— = mr[e (y) + e ±o (0)] 
k e = m[ei(y) + r(e ±1 (0) + e ±1 (0))] (18) 

Here we have introduced the vectors eo = dei/dy and e±o = dej_i/d(f> (note that all the 
occurring vectors are orthogonal). We may evidently use the quantity r (together with 
suitable boundary values) to label the wave packet for the current. That is to say, we 
may assume that there is a distribution h(r) which describes the occurrence of the different 
current lines, each with a well-defined direction r. This distribution, h(r), should be similar 
to a Gaussian. A single current line with fixed r may also be described in the transverse 
plane. The current turns around the 1-axis with the azimuthal angle and the corresponding 
field quanta are emitted transversely to the current at every emission point according to 
Eq. (p"8D . There is one reasonable restriction: the field energy emitted by the current in a 
small angular segment should not exceed the energy which should be available in the Lund 



Model string. If we use the string radius as calculated in Eq.fljJD, £ min = yTc/n, then we 
find that 

mr < Ki min ~ 0.8 GeV (19) 
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Figure 5: The break-up of a Lund string. 

It is interesting to note that these fields have similarities to those studied in connection with 
dimensional reduction in |J. 

6 Fragmentation and screwiness 

We have in the previous section described the emission of a continuous stream of colour- 
connected gluons having the property that the azimuthal angle of the stream is proportional 
to the rapidity, i.e. it is of a helical character. As previously discussed we cannot implement 
this as individual gluonic excitations of the Lund string. We will in this section instead 
describe a possible way to take the transverse properties of the continuous helix into account 
whilst keeping the major properties of the Lund fragmentation model. In order to do this 
we will begin by presenting a few relevant parts of the Lund model. This model has been 
described several times and a recent investigation can be found in [HJ . 



6.1 The Lund fragmentation process 

The following (non-normalised) probability to produce a set of hadrons has been derived 
using semi-classical arguments in |Tl"| 

n 

dP({p} f P tot ) = IlNjdpjStf} - m*)5(5>i - Ptot) exp(-bA). (20) 
i 

Here Nj are normalisation constants, A the decay area, cf. Fig.[|, and b a basic colour- 
dynamical parameter; from comparison to experimental data we know that b ~ 0.6 GeV~ 2 
if the area A is expressed in energy-momentum space quantities. 

The constant force field spanned between a colour-3 quark and a colour-3 anti-quark is 
a simple mode of the massless relativistic string. The process has been generalised into a 



situation with multigluon emission in |L2[ using the Lund interpretation that the gluons are 
internal excitations on the string field. 
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Figure 6: The longitudinal energy scale in the Lund model is (r). The figure shows the production 
of a particle with transverse mass m±. The difference in rapidity between the constituent vertices 
has to be related to the T of one of the vertices in order for the vertices to be specified. 

The area decay law in Eq. (|20f) can be implemented as an iterative process, in which the 
particles are produced in a stepwise way ordered along the positive (or negative) light-cone. 
If a set of hadrons is generated, each one takes a fraction z of the remaining light-cone 
component E + pi (or E — pi, if they are generated along the negative light-cone), with z 
given by the distribution 

f(z) = N ^~ Z ^ exp(-bm 2 Jz). (21) 

z 

The parameters N, a and b are related by normalisation, leaving two free parameters. 
The transverse mass parameter in the fragmentation function is m\ = m 2 + p]_, with the 
transverse momentum obtained as the sum of the transverse momenta stemming from the q 
and q particles generated at the neighbouring vertices, p± = k± 2 — k±_i- In the Lund model a 
qq-pair with transverse momenta ±k± is produced through a quantum mechanical tunneling 
process. It results in a Gaussian distribution for the transverse momenta 

d 2 k ± exp{-7ik 2 jK). (22) 

The whole process is implemented in the Monte Carlo program JETSET ||. 

Consider the production of a particle with transverse mass m±. Given that one vertex has 
the rapidity yi, the rapidity difference Ay is not enough to specify the position of the other 
vertex. One must also know the proper-time of the first vertex. This is shown in energy- 
momentum space in Fig|| where the first vertex is specified by T which is the squared product 
of the proper-time and k. Of course there are two solutions in this case, but one is strongly 
favoured by the area dependence in Eq. (pO[) . In the Lund model the vertices, on average, 
lie on a hyperbola given by a typical T. That is to say, the steps in rapidity in the particle 
production are related to the scale (T) as given by the model. There is a similar situation 
in the transverse momentum generation. The squared transverse momentum of a particle 
is not only given by the azimuthal angle A0 between the break-up points that generate the 
particle. The lengths of the transverse momenta of the q and the q that make up the particle 
are also needed. In the tunneling process in Eq.(^) these sizes are given by the scale k/tt. 

Thus the Lund fragmentation model provides two different energy scales; one longitudinal 
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Figure 7: A qq-pair is produced in a break-up point with azimuthal angle <f>. The figure illus- 
trates how the screwy gluon field between (p\ and fa is associated with the transverse momentum 
(—k±i, k±2) of the quarks produced at the two break-up points. This association has the property 
that the produced transverse momentum is conserved locally in each break-up point. 

to relate to the rapidity difference between vertices and one transverse to relate to their 
difference in azimuthal angle. 

6.2 A modified fragmentation process with screwiness 

The main idea in the screwiness model is that the transverse momentum of the emitted 
particles stems from the piece of screwy gluon field that is in between the two break-up 
points producing the particle. Therefore we begin by summing up the transverse momentum 
that is emitted between two points along the field line, cf. Eq.(|l8|): 



We note that the quantity mr also occurs here. We will always consider the parameter m to 
be a suitable fixed mass parameter but according to the assumed wave function for the current 
the direction r may vary between the different break-up points. To keep the presentation 
clear we will start off keeping r fixed. In the end we will present the generalisation to the 
case of a varying r. 

If we associate ±k±i with the transverse momenta of the qq-pair produced at vertex i, 
the transverse momenta of the produced particles are given by Eq. (p3|) . The corresponding 
squared transverse momentum is then 



where A0 = <pi — <pi-i. Since A0 is proportional to the rapidity difference between vertices 
Ay, it can be written as a function of the particle's light-cone fraction z 




(23) 



p\. = 2mV[l - cos(A0)] 



(24) 




(25) 
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where T is defined as in Fig|5] and with respect to the previous break-up point % — 1. Taken 
together this means that we can write the transverse momentum of a particle as a function 
of z and r 

(Ay(zY 



p 2 ± (z) = 2m 2 T 2 



1 — cos 



(26) 



As explicitly manifested in Eq.(pBl) this means that the transverse and longitudinal com- 
ponents are connected in this model. Inserting p±(z) in the Lund fragmentation function 
gives 

f{z) = ivil^l! eX p Li (ml+p 2 ± (z))\ . (27) 

In this way Eq.fl2~7|) gives the distribution of light-cone fractions for a given direction r. 

This model keeps the longitudinal properties of the ordinary Lund fragmentation model, 
but the azimuthal properties are changed. Rapidity differences are still related to (r) but 
steps in the azimuthal angle are now correlated with steps in rapidity. The azimuthal angles 
are no longer related to k/tc as given by the tunneling process, but instead to m 2 r 2 as given 
by the screwy gluon field. 

When going from one vertex to the next in the case of varying r one has to keep in mind 
that the transverse momentum produced at the first vertex has been specified by the previous 
step. In order to conserve the transverse momenta generated at each vertex we therefore 
modify the association in Eq.(^), as follows 

' Wl 1 - r *u-—L-x (28) 



-i d(f) 7;_i 

where Tj denotes the direction between break-up points i — 1 and i, cf. Fig.[7|. The transverse 
momenta of the produced particles are then given by 

p±i = k ±i - A;_Li_i = m [rie L1 (^) - Ti-ie_u((f>i-i)] , 

P±i = rn 2 [t? + rf_! - 2^-1 cos(A0)] . (29) 

The p 2 ± given by Eq. (|29|) can then be put into the fragmentation function. Varying r results 
in larger variations in the emitted transverse momenta of the particles. We have used a 
Gaussian distribution of r-directions and we have approximated m\ in Eq . (|25|) with m\ ~ 
m h + (v\min) = m h + 2m 2 o" 2 . Where m h is the hadron mass and a T denotes the width 
in the distribution of r-directions. The equations can be solved iteratively without this 
approximation, but we find that our results are unaffected by this approximation. 



7 Is screwiness observable? 

In this section we will address the question of whether introducing a correlation between y 
and of the string break-up vertices has observable consequences for the produced particles. 
There are two processes which in principle can destroy such a correlation. Firstly, there is 
the initial particle production and secondly, there are resonance decays. The initial particle 
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production spoils things because even if the vertices lie on a perfect helix the produced 
particle will usually not lie on the line between its two constituent vertices in the (y, 0)- 
plane. The particle production fluctuations are mainly in rapidity, i.e. a particle is produced 
with an azimuthal angle which roughly corresponds to the average angle of its constituent 
vertices, while its rapidity is distributed with width unity around the average of the vertices. 

To study the consequences of the screwiness model we have generated events with three 
different values (r) = 0.3,0.5 and 0.7. For each value we have tuned the parameters of 
the model to agree with the multiplicity, rapidity and transverse momentum distributions 
of default JETSET. In this way we can study the correlations introduced by the model as 
compared to the ordinary Lund string model. We have tuned m to get the default average 
of the produced particles, utilizing the fact that the product mr is the important factor. 
The parameter b has been changed from the default JETSET value to tune the multiplicity, 
and a T has been tuned to get the final charged p±_ fluctuations. Tuning with different (r) 
values results in the parameter values shown in Table [I]. We note in particular that to get 



(r) 


0.3 


0.5 


0.7 


m 


1.0 


0.71 


0.61 


b 


0.64 


0.68 


0.7 




0.2 


0.3 


0.35 



Table 1: Parameter values. The model has been tuned to the multiplicity and charged final p± 
distributions of default JETSET (b = 0.58). 



the multiplicity distributions of default JETSET only minor changes of the 6-parameter are 
needed. We also note that the restriction in Eq. ( ]T9"| ) is satisfied for all the cases since in this 
model only a fraction of the energy available in the Lund string is used to produce transverse 
momenta. 

We have generated pure qq events and the particles in the central rapidity plateau have 
been included in the analysis. The plots shown are for four units of rapidity, but the qualita- 
tive results for observable screwiness are unaffected for values as low as approximately three 
units of central rapidity. We have analysed the properties of the generated events by means 
of the screwiness measure, defined in Eq.(|13]). Here, the second sum in the measure instead 
goes over the hadrons or over the break-up vertices. The weight P e is of course unity for all 
events. 

In Fig. | the screwiness for the break-up vertices is shown. It has a clear peak for the 
different values of (r), and the a;- values for the peaks correspond to the average r values 
used. The screwiness for the initially produced particles is shown in Fig.|9]. We note that the 
peak vanishes for small values of r. A helix where the windings are separated by two units of 
rapidity corresponds to r = 1/tt. The vanishing of the signal for small r values is therefore 
in agreement with our findings for the rapidity fluctuations in the particle production. For 
comparison we have included the screwiness for the initial particles produced by default 
JETSET in Fig.[5| As expected no signal is found in this case. The screwiness is further 
diluted by resonance decays, but it is still visible for not too small r values as shown in 
Fig JO. 
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Figure 8: Screwiness for the string break-up vertices. The three curves shown are for (r) 
0.3, 0.5 and 0.7, respectively. There is a clear peak at uj ~ l/(r) in all the cases. 
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Figure 9: Screwiness for the directly produced particles. The three solid curves are for (r) = 
0.3, 0.5 and 0.7, respectively. The peak decreases as (r) is reduced. For (r) = 0.3 the peak 
has vanished due to the fluctuations in particle production. The screwiness for default JETSET 
(dashed line) has been included for comparison. 



To try to enhance the signal we have investigated how the screwiness measure depends on 
multiplicity and the transverse momentum of the particles. Selecting events with large initial 
multiplicity enhances the signal. However, analysing events with different final multiplicities 
separately does not give an enhancement of the signal. The influence of resonance decays 
on the multiplicity is too large. 

Selecting events where (p\) is large enhances the signal when decays are not included. 
This is shown for (r) = 0.3 in the left part of Fig.|ll] where events with (pj_) > 0.3 GeV 2 for 
the initial particles have been selected. As shown in the figure this event selection results 
in the signal surviving particle production even for small (r)-values. This event selection is 
also profitable when it comes to decreasing the effects of resonance decays since events with 



16 




0.075 - 



0.050 



0.025 



Figure 10: Screwiness for the final particles (7r°'s are set stable). The three curves shown are for 
(r) = 0.3, 0.5 and 0.7, respectively. For not too small (r)-values there is a peak at to ~ 1/(t). 
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Figure 11: Screwiness for (r) = 0.3. To enhance the signal events where is large have been 
selected (solid lines). We have included the corresponding curve with no event selection (dashed 
lines) to indicate the improvement. Left) Initially produced particles. {p\) > 0.3 GeV 2 . Right) 
The final particles. 7r°'s are set stable. (p 2 ± ) > 0.25 GeV 2 . 



many decay products are not likely to be selected. In the right part of Fig.[IT] we show the 
screwiness for the final state particles in events where (p 2 ± ) > 0.25 GeV 2 . The curves shown 
are for (r) = 0.3 to show that with event selection a signal can be obtained even for this 
case. Using the same event selection of course enhances the signal for larger (r)-values, but 
in those cases it was clearly visible in the total sample. 

A total of 50000 qq events have been used in the analysis, except in the event selection 



analysis in Fig.|TT] where 250000 events are analysed. To be able to observe screwiness for 
such a small (r)-value one needs to increase the number of events by a factor of about 
five compared to the larger values. Since we have only used positive (r)-values, events 
with a preferred rotational direction are generated. We could have included both rotational 
directions in the event generation which would add a signal for negative u, but reduce the 
statistics by a factor of two. 
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Figure 12: Left) Thep± (GeV) distributions for the directly produced it's (solid) and p's (dashed). 
The curves shown are for r = 0.5. Right) The p±_ (GeV) distributions for the all final pions (solid), 
7r°'s are set stable, as compared with default JETSET (dashed). 



The effects on the screwiness from hard gluons stemming from the parton cascade will be 
investigated in future work. However, since only a fairly small number of events are needed 
for the results in this paper we expect that investigations of experimental data, in which 
hard gluon activity is excluded, can be profitable. 

A specific property of our model is that (p 2 ± ) for directly produced pions is smaller than 
for heavier particles. This feature appears to be in agreement with experimental data 
on two-particle correlations |13|] . A model for correlations in p± in the string hadronization 
process with similar consequences was introduced in JT3]]. The p± for directly produced 



pions and p's are shown for the screwiness model in Fig.[L2] and the distributions are clearly 
different. In the figure we also show the p± distribution of the final pions and compare it to 
the default JETSET distribution. As seen the secondary pions wash out the differences. The 
(p±) for various flavours at the initial production level depend on the screwiness parameters, 
but the qualitative difference remains. 



8 Conclusions 



It is perhaps surprising that such an ordered structure as a helix could emerge at the end of 
the QCD cascade. However, when we consider the constraint imposed by helicity conserva- 
tion, we see that purely random configurations of gluons are disfavoured. This is because the 
exclusion region around each gluon restricts the maximum number of allowed gluons. Instead 
we see that the gluons can achieve the maximum concentration by close packing themselves 
into the form of a helix. The fragmentation of this screwy field has consequences for the 
final state particles. Although the fragmentation cannot be described in terms of gluon exci- 
tations of the Lund string, we have instead modified the Lund fragmentation scheme. If the 
winding is within reasonable limits then we expect "screwiness" to be an observable feature 
of the QCD cascade. 
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Figure 13: The multiplicity in central rapidity per unit of rapidty n and the corresponding variance 
a\ depends on the upper cut-off in the cascade &_i_ m ax (GeV) as shown. Default JETSET has been 
used for fragmentation and cr^/n does not start to increase until k± max is roughly 1.6 GeV 
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A Problems with fragmenting soft gluons 

In section |] we claimed that the helix colour-field cannot be implemented as an excited 
string, since gluons softer than k±o = 1.6 GeV cannot be considered as excitations of the 
string. 

To illustrate the problems with fragmentation of soft gluons we have investigated JETSET 
fragmentation of parton configurations with soft gluons emitted according to the Dipole 
Cascade Model as implemented in the ARIADNE Monte Carlo flSfl . The allowed k± range 



for emissions from the colour dipoles is normally between an upper value, given by phase- 
space limits, and a lower infra-red cut-off, k^ c . We have instead used a small maximum 
allowed k± value (denoted fc_i_ max ) to restrict the hardness of the emitted gluons. This soft 
cascade has been applied to qq-dipoles oriented along the z-axis. The soft gluons have a 
negligible impact on the event topology and for our purposes it therefore makes sense to 
define rapidity with respect to the z-axis. We have analysed the resulting hadrons in the 
central rapidity plateau of the events. To emphasize the features of fragmentation of soft 
gluons we have not included resonance decays in our analysis. 

In Fig.|T3| we show how the average and the squared width of the central multiplicity 
distribution depend on /c_i_ max - The effect of the soft gluons is an increase of the average 
multiplicity while the multiplicity fluctuations remain constant or even decrease until fc± ma x 
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is above kxo- The (p±) with respect to the z-axis of the hadrons only increases from 0.46 
GeV for a flat string with no gluon excitations to 0.56 GeV for k± max = 3 GeV. Changing the 
generated (p±) by such a small factor has a minor effect (~ 5%) on the average multiplicity 
in pure qq events whilst adding the soft gluons increases the average multiplicity by roughly 
40%, as shown in the figure. As mentioned in section^, we find that the soft gluons essentially 
only increase the hadron multiplicity. The number of gluons per rapidity unit varies from 
0.25 for fcj_ max = 1 GeV to 0.7 for k± max = 5 GeV. The situation is even worse in the case of 
the helix field where the expected number of soft gluons per unit of rapidity is significantly 
larger. We conclude that gluons softer than k^ cannot be implemented inside the Lund 
Model as individual gluonic excitations of the string. 

We will end this appendix with an interpretation of the Lund fragmentation model, which 
provides us with the possibility to relate kj_o to the 6-parameter in the model. The result 
in Eq.(|20|) (although derived semi-classically) can be interpreted quantum-mechanically by 
a comparison to Fermi's Golden Rule. It equals the final state phase space times the square 
of a transition matrix element \Ai\ 2 = exp(— bA). There are two such quantum-mechanical 
processes, Schwinger tunneling and the Wilson loop integrals, which can be used in this 
connection (and they result in very similar interpretations of the parameters). For the 
Schwinger tunneling case we note that if a constant (k) force field is spanned across the 
longitudinal region X during the time T with a transverse size A± then the persistence 
probability of the vacuum (i.e. the probability that the vacuum should not decay by the 
production of new quanta) is [|16| 

\M\ 2 = exp{-^ 2 XTA L U) . (30) 

Here the number II only depends upon the properties of the quanta coupled to the field; for 
two massless spin 1/2 flavours it is II = 1/127T. Comparing the result in Eq.(|3~0"l) to Eq. ([20|) 
we find that the parameter b = A±/2Att (taking into account that the Lund model area is 
counted in lightcone units). From Eq.(|l4]) we obtain the minimum transverse size of the field 
from which it then follows that the 6-parameter in the Lund model must be b > 7t/24k ~ 0.6 
GeV -2 . This is evidently just in accordance with the phenomenological findings in the Lund 
model for the parameter b. Further, considering the distribution in Eq. (f22|) for the transverse 



momentum of a produced qq-pair breaking the string we recognise the quantity in the 
exponential fall-off. We may conclude that there is a wave-function for the Lund string 
in transverse space with just the right transverse size to allow the "ordinary" transverse 
fluctuations in momenta. 
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